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HIGHLIGHTS 


• A semi-analytical method for the HPFE parameter design is proposed. 

• The HFPE piston oscillation is a hard self-oscillation. 

• The piston oscillation has an unsymmetrical oscillation characteristic. 

• Stable running of a single piston HFPE prototype is realised. 
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The hydraulic free-piston engine is a hybrid of a hydraulic plunger pump and a free-piston engine. This 
paper presents a semi-analytical method for the hydraulic free-piston engine parameter design. An oscil¬ 
lation function described the piston oscillation were established. The piston oscillation characteristics 
were analysed. The parameter sensitivity and the limit cycle stability of the piston oscillation were inves¬ 
tigated. The results show that the single piston free-piston engine system is a hard self-oscillation with 
variable damping and stiffness. External energy is required to start the piston oscillation. The piston oscil¬ 
lation has an unsymmetrical oscillation characteristic. The dead centres of the piston motion are deter¬ 
mined by the energy balance. The energy supply and the energy consumption are determined by the 
generalised damping in the piston oscillation function. The piston oscillation frequency should satisfy 
the minimum compression ratio requirement for the compression ignition. In the steady operation, the 
piston motion is a self-excitation oscillation. The trajectory stability is mainly determined by the fuel 
combustion energy supply. The scavenging affected by the BDC should be paid special attention to ensure 
the trajectory stability. 

© 2014 Elsevier Ltd. Ah rights reserved. 


1. Introduction 

The free-piston concept was firstly proposed by Pescara in 
1920s [1] and developed fast from 1930s to 1960s [2,3]. However, 
for many reasons, the free-piston engine (FPE) was replaced by ad¬ 
vanced crank internal-combustion engines and practically aban¬ 
doned since the 1960s [3 . Potential advantages of the FPE, such 
as optimised combustion through variable compression ratio, and 
lower frictional losses due to the simple structure with few moving 
parts, makes the FPE concept receive more interest in recent years 
[4-7]. 

Compared to the conventional crank internal-combustion en¬ 
gine, the FPE has two typical features [2,3]: 

(a) The piston motion is not restricted by a rotating crankshaft, 
as known from crank engines. 
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(b) The power is taken out by media, not by gear or shaft. 

The FPE concept can be divided into three categories according 
to the number or position of the piston. With different media of 
loads, the FPE can be converted into different power engines. Most 
of them are designed as power sources for automotive applications 
[3 . The General Motors XP-500 completed in 1957 was the first 
free-piston automobile [3]. The XP-500 was powered by the Gen¬ 
eral Motors Research 4-4 ‘Hyprex’ engine. The General Motors Re¬ 
search 4-4 ‘Hyprex’ was a dual opposed piston, diesel free-piston 
engine. It had a power output of around 185 kW [8 . Because of 
technical shortcomings and more promising research in other 
areas of alternative propulsion, General Motors halted further 
development on the free-piston automobile within three years of 
the XP-500’s completion 9]. Ford, like General Motors, saw a po¬ 
tential in the free-piston gas generator for automotive applications 
in the mid-1950s [3 . Frey and co-researchers proposed an analyt¬ 
ical model for the design of a free-piston gas generator and de¬ 
scribed the development of a 112 kW prototype [10. The 
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Nomenclature 



a 

amplitude 

Pgc 

gas pressure brought by the fuel combustion 

A, 

function of the amplitude 

Pgv 

gas pressure in the pure compression process 

b 

decentration 

Pw 

hydraulic working pressure 

c 

generalised damping 

Para 

design parameter value 

Cequi 

equivalence damping 

Pal CZ rated 

design parameter rated value 

Co 

damping coefficient 

S 

piston cross-section area 

D 

cylinder bore 

Si 

cross-section area of the check chamber 

E 

^cyc 

energy released by fuel combustion in one cycle 

S 2 

cross-section area of the pump chamber 

/ 

piston oscillation frequency 

S 3 

cross-section area of the compression chamber 

fc 

simplified combustion gas force 

Vo 

initial value of the in-cylinder gas volume 

F C h 

check chamber force 

y 

state variable 

Fco 

compression chamber force 

X 

piston displacement 

Fg 

in-cylinder gas force 

X* 

displacement in the calculation coordinate system 

Fgc 

force brought by the fuel combustion 

x\ 

function of the amplitude 

fgCmax 

maximum value of force brought by the fuel 

*CSBDC 

BDC in the compression stroke 


combustion 

*ESBDC 

BDC in the expansion stroke 

F g v 

force determined by the changing of gas volume 



F pu 

pump chamber force 

Greek symbols 

F r 

friction force 

y 

in-cylinder gas adiabatic exponent 

F r 0 

value of friction force 

3 

active range of the combustion force 

I< 

generalised nonlinear stiffness 

s 

compression ratio 

L 

distance between the piston balance position and the 

Pi 

indicated thermal efficiency 


cylinder head plate 

(p 

phase 

m 

moving piston mass 

<P> 0 and functions of the amplitude and the phase 

P 

in-cylinder gas pressure 

CO 

change rate 

Po 

initial value of the in-cylinder gas pressure 



Pc 

hydraulic compression pressure 




favourable torque-speed characteristic of the prototype led Ford 
engineers to conclude that the Ford model 519 free-piston power 
plant was ideally suited for installation in the Ford Typhoon tractor 
[3]. A hydraulic free-piston engine (HFPE) has been proposed as a 
low-cost and high-efficiency power unit for a fork-lift truck by In¬ 
nas [4]. The FIFPE is a combination of a free-piston engine and a 
hydraulic pump. The HFPE gets higher efficiency compared to the 
conventional assemble of a crank internal-combustion engine 
and a hydraulic plunger pump [11]. 

The free-piston engine characteristic is often investigated by 
discrete numerical methods [12-22,7 . Zero-dimensional model¬ 
ling of free-piston engines has been proposed [12-19 . Some have 
investigated free-piston engines using multidimensional simula¬ 
tion models [20-22,7]. The numerical methods could achieve the 
instantaneous displacement, velocity and acceleration of the en¬ 
gine piston with a high degree of precision. However, the numeri¬ 
cal result only contains the discrete numerical solution and the 
panorama of the system solution is still not completely offered. 
To achieve high precision numerical result, the numerical model 
is always complicated and the time-consuming for the model cal¬ 
culation is long. The numerical result is inconvenient for analysing 
the piston oscillation and the contributing factors. The result also 
cannot fully satisfy the requirement for the engine parameter con¬ 
trol system design. A semi-analytical model has been developed for 
the free-piston Stirling engine analysis [23]. 

This paper studies the design method of a single piston hydrau¬ 
lic free-piston engine (SPHFPE) as the power source of a hydraulic 
hybrid vehicle. Different from the discrete numerical method used 
in the engine parameter design, a semi-analytical method was 
used. The forces applied on the piston were analysed and the pis¬ 
ton oscillation semi-analytical model was established. Extensive 
results of the oscillation model were presented, giving insight into 
piston oscillation characteristics. The transient process between 
different stable piston trajectories was analysed. The results 


indicate that the steady state of the designed free-piston engine 
is a hard self-excited vibration with variable damping and stiffness. 
The stable piston trajectories only appear in a special piston stroke 
range, which means that the engine is operated in energy balance 
conditions. The aim is to offer the piston oscillation panorama by 
the semi-analytical method. 

2. Configuration of the SPHFPE 

The configuration of the SPHFPE is shown in Fig. 1, firstly pro¬ 
posed by Innas BV [4 . The SPHFPE is a two-stroke diesel engine 
with uniflow scavenging and direct fuel injection. The exhaust 
valves and the scavenging pump are driven by the high-pressure 


Fig. 1 . SPHFPE configuration. CD Combustion cylinder, © check chamber, CD pump 
chamber, © compression chamber, © reset valve, © piston displacement trans¬ 
ducer, CZ) frequency control valve, ® compression accumulator, ® high-pressure 
rail, © low-pressure rail. 
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^Compression accumulator 
Frequency control valve 
/ Fuel injector 


Fig. 2. SPHFPE prototype. 


rail. For the injection system, a hydraulic-electronic unit injector 
from Caterpillar is used, as shown in Fig. 2. 

The hydraulic pump of the SPFIFPE consists of the check cham¬ 
ber, the pump chamber and the compression chamber. The check 
chamber restricts the rebound of the piston around bottom dead 
centre (BDC) after the expansion stroke. The pump chamber deliv¬ 
ers one part of the high pressure oil in the expansion stroke. The 
other is delivered by the check chamber in the compression stroke. 
The compression chamber completes the compression stroke by 
the hydraulic power stored in the compression accumulator. The 
frequency control valve is used to control the working frequency 
of the SPFIFPE. The design parameters of the SPHFPE prototype are: 

• Piston diameter: 98.5 mm. 

• Stroke: 114-125 mm. 

• Mass of the moving piston: 4.2 kg. 

• Pressure of the high-pressure rail: 20-30 MPa. 

• Pressure of the low-pressure rail: 1.0-2.0 MPa. 

• Hydraulic compression pressure: 15-16 MPa. 

• Maximum output power: 15 kW. 

3. Forces on the piston 

The forces on the piston include in-cylinder gas force, hydraulic 
forces and friction forces, as shown in Fig. 3. The in-cylinder gas 
force F g is brought by the in-cylinder gas pressure p. The force 
can be divided into two components. The force F gV is determined 
by the changing of gas volume. The force F gC is brought by the fuel 
combustion. The relation between the two gas forces is 

= -Fgv + F g c (1) 

The system coordinate was chosen with the origin at the cylin¬ 
der head plate and the positive direction points to BDC. A dynamic 
model of the piston motion was developed based on Newton’s sec¬ 
ond law: 

mx = F g + F ch - F pu - F co - F r sign(x) - c 0 x (2) 

where m is moving piston mass, x is piston displacement, F r sign(x) 
is friction force, c G is constant damping coefficient. 

Substituting (1) into (2), we can get 

mx = F gV + F gC + F ch - F pu - F co - F r sign(x) - c 0 x (3) 

Eq. (3) can be rewritten as: 

mx + c(x)x + I((x 7 x)x = 0 (4) 



where 

c ( . )= Fp U+ F rS ign(x)-F 8C + Co (5) 

X 

K(x,x)= Fch+Fco ~ FgV (6) 

X 

Eqs. (3)-(5) describe a single degree freedom nonlinear vibra¬ 
tion system. The c(x) represents a generalised damping, deter¬ 
mined by the friction force, the pump chamber force, the 
combustion gas force and the piston velocity. The JC(x,x) represents 
a generalised nonlinear stiffness, determined by the check cham¬ 
ber force, the compression chamber force, the volume changing 
gas force and the piston displacement. The generalised nonlinear 
stiffness does not consume or supply energy in one engine cycle. 

The signs of the generalised damping items varying with the 
piston displacement are shown in Fig. 4. The plus sign means 
absorbing the piston kinetic energy, while the minus sign means 
releasing. The piston displacement is asymmetric according to 
top dead centre (TDC) [24]. The friction force part of the general¬ 
ised damping is always to be positive in engine cycles, which 
means that it consumes the piston kinetic energy all the time. 
The sign of the pump chamber force changes to be the plus sign 
when the piston crosses TDC. This indicates that the piston kinetic 
energy is absorbed by the pump chamber force in the expansion 
stroke. The combustion gas force appears only after the injection 
of fuel. The sign of the combustion gas force part changes to be 
negative. It means that the energy is supplied by the fuel combus¬ 
tion and the piston kinetic energy is increasing. Further, if the com¬ 
bustion gas force appears in the compression stroke, it results in a 
decrease of the compression ratio. 

To make the piston moving continually, the energy consumed 
and supplied by the generalised damping should be equal and en¬ 
ough. Based on the piston motion characteristics, the piston oscil¬ 
lation is a self-oscillation. When the energy supplied by the 
combustion gas force increases, the piston stroke also rises cycle 
by cycle. With the increase of the piston stroke, the energy con¬ 
sumed by the friction force and the pump chamber force also in¬ 
creases until the energy balance is achieved again. When the 
releasing energy decreases, the piston stroke falls until the energy 
balance is achieved in the same way. Namely, if the energy sup¬ 
plied by the combustion gas force is provided stably cycle by cycle, 
the piston reciprocating motion continues. If the supplied energy 
keeps being constant, the piston stroke is also stable and the piston 
oscillation becomes steady. The piston steady oscillation is a self- 
excited vibration. The steady phase trajectory of the piston is a lim¬ 
it cycle, which is a closed trajectory in the phase space. For the 
nonlinear vibration of the piston, the limit cycle means that at least 
one other trajectory of the piston spirals into it as time approaches 
infinity. The piston stroke becomes stable and the energy balance 
is achieved at last. 


X 

Item 
F> ign(x) 

Compression stroke Expansion stroke 


+ + + + + + + + + + 

X 

F PU 


+ + + + 

X 

5c 


X 




Time t 


Fig. 3. Forces applied to the piston. 


Fig. 4. Signs of the generalised damping items. 
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4. Piston oscillation model 

The piston oscillation characteristics are analysed by the phase 
plane method: 

x 

c(x)x + K(x)x (7) 

m 

where y is state variable. 

For further understanding of the system, the analytical method 
is used to get the piston limit cycle. The check chamber force and 
the compression chamber force are assumed to be constant. The 
pump chamber force is only applied on the piston in the expansion 
stroke. Eqs. (8)-(10) give the expressions of the forces: 


F ch — P w Sl 

(8) 

Fco = Pc$ 3 

(9) 

p w S 2 sign(x)+p w S 2 

F pu = 2 

(10) 


where p w is hydraulic working pressure, Si is cross-section area of 
the check chamber, p c is hydraulic compression pressure, S 3 is 
cross-section area of the compression chamber, S 2 is cross-section 
area of the pump chamber. 

The volume changing gas force is calculated by 

Fgv=Po(y)'s (11) 

where p 0 is initial value of the in-cylinder gas pressure, V 0 is initial 
value of the in-cylinder gas volume, S is piston cross-section area, y 
is in-cylinder gas adiabatic exponent. 

As shown by Fig. 5(a), p gV is in-cylinder gas pressure in the pure 
compression process and p gC is gas pressure brought by the fuel 
combustion. In order to simplify the complex fuel combustion pro¬ 
cess, the force applied by fuel combustion is replaced by the force 
shown by Fig. 5(b). The force is also assumed to actuate only in the 
expansion stroke. The simplified combustion force is easily de¬ 
scribed by a simple mathematical expression and very useful for 




BDC 

TDC 


Time t 


Fig. 5. The simplified combustion gas force, (a) Gas pressure components, (b) 
Simplified method. 


raising the model simulation speed. It seems that proper precision 
is also achieved, as shown in Fig. 5. The simplified combustion gas 
force can be expressed as: 

F gC =/ c (x,x) (12) 

Based on the energy conservation law, the force maximum 
value Fgcmax should satisfy: 

(13) 

where F cyc is energy released by fuel combustion in one cycle, ^ is 
indicated thermal efficiency, <5 is action length of the force Fgo 
The friction force is 

FfSign(x) = Frt,sign(x) (14) 

where F r0 is value of the friction force. The value is a constant. 

To get the limit cycles, a system of coordinates was chosen 
having the origin at the piston moving balance position and the 
positive direction pointing to TDC. Then Eq. (3) can be transformed 
to: 


mx* + 


7ID Po 


X* 


1+T| J ) + !frtli(|) i + w+')0 + 2 ) (|' 1 


+/c(x*,x*)+p w ( 5 i~y) + 


(2F r + p w S 2 )sign(x*) 


-p c S 3 =0 (15) 


where x* is displacement in the calculation coordinate system, D is 
cylinder bore, L is distance between the piston balance position and 
the cylinder head plate. 

With a simplified method, Eq. (15) is represented as 


*‘+g(x*)=F/(x*,x‘) 


(16) 


where 


g(x*) 


PoHf)' 

ml 


x* + 


(1 + y)x* 

21 


(17) 


/(x*,x*) = - 


f c (x‘,x')+p 0 n 


D‘ 




2 , F r Jsign(x*)+p 


W 


Si 1 -p c S 3 


(18) 


The first order approximate solution of Eq. (16) and the approxi¬ 
mate solution is given by 

( x* = a cos (p + b + (a) 

a = hMa) (19) 

<P = <P 0 (a,(p)+^( a ^) 


where a is amplitude, b is decentration, cp is phase, x\ and A\ are 
functions of the amplitude, & 0 and are functions of the 
amplitude and the phase. 


5. Results and discussion 

5.2. Model validation 

The measurement system of the prototype had pressure mea¬ 
surements on the combustion cylinder, the hydraulic chambers 
and in the common pressure rails, as shown in Fig. 6. The used 
measuring sensors are presented in Table 1. The piston displace¬ 
ment has been measured directly by a laser displacement sensor, 
since there is no mechanisms coupled to the SPHFPE piston assem¬ 
ble. The measured data synchronisation of the combustion cylinder 
pressure and the piston displacement is important and the data 
acquisition system PXI-6123 has been used. Further, the maximum 
piston acceleration appears at the same time with the highest com- 
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Displacement 

Pressure 



-(Charge amplifier)- 



DAQ system 


Fig. 6. Schematic diagram of experimental test rig. 


Table 1 

Used measuring sensors. 


Measuring parameter 

Sensors 

Combustion cylinder pressure 
Hydraulic chamber pressures 
Piston displacement 

Exhaust valve stroke 

Engine body axis displacement 

Kistler 6055BB, Charge Amplifier 501 IB 
HYDAC HDA-3844 

Laser displacement sensor ZLDS100/HS 
Laser displacement sensor LD1607 
LVDT-EK1000 



Piston displacement x / m 

Fig. 7. Measured and calculated piston motion under given parameters. 


bustion cylinder pressure should be mentioned. The information 
was collected by National Instruments data acquisition system 
PXI-8196 and PXI-6123. The response of the data acquisition sys¬ 
tem is 50 kHz and the frequency has been applied to all the mea¬ 
sured parameters individually. Further, a DAQ-board 6259 from 
National Instruments is also used to collect both the control signals 
and the sensor information. The response of the data acquisition 
system is 5 kHz individually. 

The operation parameter of the SPHFPE prototype is presented 
in able 2. The calculated piston oscillation was compared to data 
from the SPHFPE prototype, as shown in Fig. 7. The comparison 
was done with the aim of verifying that the model produces real¬ 
istic results, and that it is able to predict real trends for varying en¬ 
gine design parameters. 

5.2. Transient process 

Before the SPHFPE piston oscillation takes place, initial oscilla¬ 
tion energy is needed. The initial oscillation energy provides the 
gas compression energy of the start engine cycle and is stored in 
the compression accumulator. Without the initial oscillation en¬ 
ergy, the SPHFPE piston oscillation never starts by itself. It seems 
that the engine piston dynamic system is a hard self-excited oscil¬ 
lation system. 

The phase locus of different start compression positions is 
shown in Fig. 8. The piston oscillation is an unsymmetrical oscilla¬ 
tion. It seems that the start oscillation position affects the phase lo¬ 
cus of the engine little. However, the compression ratio of the start 
engine cycle should be higher enough since the SPHFPE is a diesel 
engine. The initial oscillation energy stored in the compression 
accumulator is transformed into the gas internal energy in the start 
engine cycle. It is seen that the piston oscillation becomes stable 
with the constant fuel combustion energy and the constant 
hydraulic working pressure. For the SPHEPE, the energy released 
by the fuel combustion is important for the piston oscillation sta¬ 


bility. The released energy determines the expansion stroke length 
since the hydraulic working pressure changes little. The variable 
expansion stroke changes the BDC. The flow area of the intake port 
is affected by the BDC changes. Thus, the scavenging effect is 
mostly determined by the BDC under the constant scavenging 
gas pressure. The released energy needs to make the piston move 
back to the BDC. Then, the intake port is opened widely enough. 
The piston oscillation limit cycle is determined by the energy bal¬ 
ance between the fuel combustion energy and the hydraulic load 
energy. It seems that a limit cycle appears in any piston stroke if 
the stable fuel combustion energy is supplied. 

In the steady operation, the piston oscillation is a periodic oscil¬ 
lation. An equivalence damping can be used instead of the general¬ 
ised damping in the steady operation. The calculation of the 
equivalence damping is based on the equality of the energy con¬ 
sumed by the generalised damping. The equivalence damping is 
calculated by 

C X ESBDC T X ESBDC 

/ c(x)dx = / c equi (x)dx (20) 

7xcsbdc 7x C sbdc 

where x C sbdc is BDC in the compression stroke, x ES bdc is BDC in the 
expansion stroke, c equ i(x) is equivalence damping. 

Then, Eq. (4) can be transformed to 

mx + X(x,x)x = 0 (21) 

Base on Eq. (21), an amplitude frequency response characteris¬ 
tic of the equivalence oscillation system is given. The compression 
ratio s which is especially important for the compression ignition 
engine is confirmed. The relationship between the compression ra¬ 
tio and the piston oscillation frequency/is presented based on Eq. 
(21), as shown in Fig. 9. The designed piston moving mass of the 
SPHFPE is about 4.2 kg. It is seen that the piston oscillation fre¬ 
quency has to satisfy the minimum compression ratio requirement 
for the compression ignition. The compression ratio increases with 
the increase of the piston oscillation frequency. Because the output 
power is in direct proportion to the piston oscillation frequency, an 


Table 2 

Used control parameters. 


Parameter 

Value 

Parameter 

Value 

Output pressure 

12.8 MPa 

Low-pressure rail 

1.2 MPa 

Compression pressure 

14.0 MPa 

Injector driving pressure 

14.0 MPa 

Injection timing (before TDC) 

23 mm 

Valve timing (after TDC) 

58 mm 

Exhaust valve driving pressure 

14.0 MPa 

Exhaust valve signal pulse width 

10 ms 
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Piston displacement x / m 


Fig. 8. The phase locus under different piston strokes, (a) A phase locus with 



Piston displacement x / m 


a decreasing piston stroke, (b) A phase locus with an increasing piston stroke. 




Oscillation frequency// Hz Piston displacement x / m 

Fig. 9. Compression ratio vs. oscillation frequency of different piston masses. Fig. 10. The phase locus of different piston moving masses. 


effective method is required to change the output power of the 
SPHFPE. To get a lower piston oscillation frequency, a pulse pause 
modulation control method has been proposed [4 . The pulse 
pause modulation control does not change the piston oscillation 
characteristics of one cycle. Thus, the compression ratio require¬ 
ment and the controllable output power are satisfied at the same 
time. 


5.3. Parameter sensitivity 

The piston oscillation characteristics change with different en¬ 
gine design parameters. To indicate the changing characteristics 
of the piston motion more clearly, the relation between the piston 
displacements x and x* is given by 


f x = I - x* 
\ x = -x* 



Fig. 10 shows the piston oscillation trajectory at the different 
piston moving mass. It seems that the change of the piston moving 
mass have little effect on the engine dead centre. The different pis¬ 
ton moving mass only leads to changes in both the piston velocity 
and the piston oscillation frequency. The engine gets a higher pis¬ 
ton oscillation frequency with a smaller piston moving mass. How¬ 
ever, the higher piston oscillation frequency demands a much 
faster response of the engine control system. The exhaust valve 
needs a fast control unit, which has been designed for the SPHFPE 


[25]. The spool of the check valve on the pump chamber also need 
to move fast to reach a higher volumetric efficiency. Some methods 
have been proposed to satisfy the fast response requirement. 

Fig. 11 proposes the piston oscillation trajectory and the dead 
centre change at the different hydraulic working pressure. The 
change rate of a design parameter is calculated by 


co = 100 x 


Para - Para ratec i 
Para rated 



where Para is design parameter value, Para ratec i is design parameter 
rated value, co is change rate. It is seen that the change of the 
hydraulic working pressure has a direct effect on the piston stroke 
because of the energy balance in the engine operation. Both TDC 
and BDC have shown a high sensitivity to the hydraulic working 
pressure. With the increase of the hydraulic working pressure, the 
TDC position increases and the compression ratio becomes smaller. 
The change rate of the compression ratio is about two times larger 
than the hydraulic working pressure change rate. The increase of 
the hydraulic working pressure may lead to misfire since the in¬ 
crease of the hydraulic working pressure has a great effect on the 
decrease of the compression ratio. The increase of the hydraulic 
working pressure also causes the reduction of the BDC. The piston 
stroke also decreases. It seems that a stable hydraulic working pres¬ 
sure is favourable for the engine operation. The minimum compres¬ 
sion ratio of the engine should be guaranteed in the engine design 
when the maximum hydraulic working pressure is applied. Further, 
the measured hydraulic working pressure needs to be set as a load 
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Piston displacement x / m 
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Fig. 11. The phase locus of different hydraulic working pressures, (a) Phase locus variety, (b) Dead centre changes, (c) Values of the dead centres. 
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Fig. 12. The phase locus of different hydraulic compression pressures, (a) Phase locus variety, (b) Dead centre changes, (c) Values of the dead centres. 


feedback for the engine control unit design, since the load of the 
engine has a linear relationship with the pressure. 

Fig. 12 shows the piston oscillation trajectory of different 
hydraulic compression pressures. It can be seen that the hydraulic 


compression pressure has a much larger effect on the TDC than the 
BDC. The piston is free to move linearly and the compression stroke 
is realised by the hydraulic energy provided by the hydraulic accu¬ 
mulator. The hydraulic compression pressure in the hydraulic 
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accumulator is an important parameter for the compression ratio 
control. An approximately linear relationship between the hydrau¬ 
lic compression pressure and the TDC is shown. Thus, the hydraulic 
compression pressure largely determines the compression ratio. 
The compression ratio can be adjusted freely by controlling the 
hydraulic compression pressure. The engine can also get a better 
compression ignition performance by increasing the hydraulic 
compression pressure, especially in the start engine cycle. Further, 
the piston stroke length changes little when the hydraulic com¬ 
pression pressure changes. To get an accurate compression ratio 
control, the changing of the hydraulic compression pressure needs 
an additional hydraulic control system. However, the minimum 
compression ratio of the engine should be achieved in the com¬ 
pression stroke with the minimum hydraulic compression pres¬ 
sure. Achten and co-researchers [4] have proposed a method to 
change the hydraulic compression pressure in the compression 
accumulator by means of electrohydraulic valves. The hydraulic 
compression pressure represents the compression energy, because 
the gas precharge pressure of the compression accumulator is con¬ 
stant. The amount of compression energy and subsequently the 
compression ratio can be controlled by the method. The changes 
of the hydraulic working pressure may change the BDC of the next 
cycle. It means the compression energy changes if the hydraulic 
compression pressure is constant. The compression ratio becomes 
variable, which is unfavourable for the engine stable operation. For 
the variations in the hydraulic working pressure, the hydraulic 
compression pressure of the next cycle should have the same 
change trend with the hydraulic working pressure. It is aimed to 
achieve a stable compression ratio under a variable piston stroke. 

Fig. 13 presents the piston oscillation trajectory changes at the 
different fuel combustion energy. It can be seen that the expansion 
stroke length is deeply affected by the fuel combustion energy. The 
stroke length is a result of the energy balance between the fuel 
combustion energy and the hydraulic load energy in the expansion 
stroke. The smaller fuel combustion energy causes a smaller BDC. 


D 



Fig. 14. Special piston positions. A is the position of the cylinder head plate. B is 
determined by the minimum compression ratio. C is determined by the smallest 
flow area of the scavenging port. The length between C and D is determined by the 
piston rebound. E is determined by the hydraulic pump end closure. 


The TDC becomes larger and the compression ratio becomes smal¬ 
ler. The trend is bad for the stability of the engine operation. The 
accurate description of the combustion heat releasing is still a 
big challenge in the free-piston engine researches. Further, the 
combustion heat releasing is also rapid and the short combustion 
duration decreases the rated engine power [24]. The smaller rated 
engine power has a negative effect on the FPE application field 
development. Higher fuel mass flow for a faster combustion of 
one cycle is needed to increase the rated engine power. In the 
SPHFPE control system, a fuel injection control strategy is used 
with the piston displacement feedback. When the BDC becomes 
too small, a little larger fuel mass is injected. The control is to make 
the piston go back to the rated BDC reliably. 


5.4. Limit cycle stability 

The piston oscillation limit cycle appears when there is suffi¬ 
cient combustion energy released. To realise the fuel combustion, 
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Fig. 13. The phase locus of different combustion energies, (a) Phase locus variety, (b) Dead centre changes, (c) Values of the dead centres. 
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Fig. 15. Measured piston oscillation characteristics, (a) Stable piston oscillation, (b) Unstable piston oscillation and misfire. 


an enough compression ratio is important for a compression igni¬ 
tion engine. The compression ratio requirement for the SPHFPE 
compression ignition is not satisfied in most small piston displace¬ 
ments. Too small piston displacement also affects the SPHFPE oper¬ 
ation, especially the scavenging effects. Guaranteeing the 
scavenging effect can improve the engine operation performance. 
To achieve a better scavenging effect, the sufficient flow area of 
the scavenging port is required. Thus, the minimum BDC should 
be confirmed, as shown in Figs. 14 and 15. The actual BDC should 
be controlled to be big enough to ensure suitable flow areas of the 
scavenging ports. Further, the scavenging effect is determined not 
only by the flow area of the scavenging port, but also by the scav¬ 
enging pressure. The higher scavenging pressure increases the 
scavenging performance. However, the higher hydraulic compres¬ 
sion pressure is also needed to achieve the same compression ratio 
under the higher scavenging pressure. For the SPHFPE, the maxi¬ 
mum BDC is restricted by the hydraulic pump end closure. Further, 
if the piston moves far outside the rated BDC, it leads to a mechan¬ 
ical contact between the piston and the end closure head. The 
mechanical contact damages the mechanical structure of the 
SPHFPE. 

In the practical application, the SPHFPE often works under its 
self-exacted frequency since the output power is not always the 
largest. Therefore, the piston oscillation frequency needs to be con¬ 
trolled. The control is realised by stopping the piston around BDC 
for a moment. The piston oscillation characteristics will not be 
changed under the control, because the piston oscillation is a hard 
self-oscillation. The control parameters are included in the gener¬ 
alised damping and the generalised nonlinear stiffness. The scav¬ 
enging pressure, the injection timing, the injected fuel mass, the 
hydraulic working pressure and the hydraulic compression pres¬ 
sure are included in the control parameters. The oscillation results 
propose the right direction to the appropriate parameter set of the 
SPHFPE control system as soon as the operation condition of the 
engine changes. Other contributing factors such as the hydraulic 
leakage and the hydraulic valve dynamic should also be considered 
in the engine operation. The effects of these factors could be 
achieved by both the full-cycle model simulation and the proto¬ 
type test. 


6. Conclusions and future work 

The semi-analytical method for the SPHFPE parameter design 
was proposed. The results indicate that the piston oscillation has 
an unsymmetrical oscillation characteristic. The oscillation is a 
hard self-oscillation with variable damping and stiffness. The 
external energy is required to start the piston oscillation. The sta¬ 
bility of the piston oscillation is determined by the energy balance 
between the fuel combustion energy and the hydraulic load. The 


compression ratio increases with the higher piston oscillation fre¬ 
quency. The piston oscillation frequency should satisfy the mini¬ 
mum compression ratio requirement for the compression 
ignition. In the steady operation, the dead centres of stable piston 
trajectories are mainly determined by the hydraulic working pres¬ 
sure, the hydraulic compression pressure and the fuel combustion. 
With the increase of the hydraulic working pressure, the piston 
stroke decreases and the compression ratio becomes smaller. An 
approximately linear relationship appears between the hydraulic 
compression pressure and the TDC. The expansion stroke is deeply 
affected by the fuel combustion energy. The scavenging affected by 
the BDC should be paid special attention to ensure the stable pis¬ 
ton oscillation. 

The results show the right direction to the appropriate param¬ 
eter set of the SPHFPE control system as soon as the operation con¬ 
dition changes. High reliability and fast response of engine control 
system are significant for the SPHFPE piston motion control. Much 
development work remains before the SPHFPE could be used for 
the realistic automotive application. Work on these topics is cur¬ 
rently underway in the National Key Laboratory of Vehicular 
Transmission at the Beijing Institute of Technology. 
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